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Several experiments reveal the inhomogeneous character of the superconducting state that oc¬ 
curs when the carrier density of the two-dimensional electron gas formed at the LaXOa/SrTiOa 
(X=A1 or Ti) interface is tuned above a threshold value by means of gating. Re-analyzing previous 
measurements, that highlight the presence of two kinds of carriers, with low and high mobility, we 
shall provide a description of multi-carrier magneto-transport in an inhomogeneous two-dimensional 
electron gas, gaining insight into the properties of the physics of the systems under investigation. 

We shall then show that the measured resistance, superfluid density, and tunneling spectra result 
from the percolative connection of superconducting “puddles” with randomly distributed critical 
temperatures, embedded in a weakly localizing metallic matrix. We shall also show that this sce¬ 
nario is consistent with the characteristics of the superconductor-to-metal transition driven by a 
magnetic field. A multi-carrier description of the superconducting state, within a weak-coupling 
BCS-like model, will be finally discussed. 

PACS numbers: 71.70.Ej,73.20.-r,73.43.Nq,74.81.-g 


I. INTRODUCTION 


After a two-dimensional electron gas (2DEG) was de¬ 
tected at the interface between two insulating oxides?® 
an increasingly intense theoretical and experimental in¬ 
vestigation has been devoted to these systems. The 
properties of this 2DEG are intriguing within several re¬ 
spects: it can be made superconducting when the carrier 
density is tuned by means of gate voltage (see Fig.[l]) , 
both in LaAlC> 3 /SrTi 03 (henceforth, LAO/STO) 12 and 
LaTi0 3 /SrTi0 3 (henceforth, LTO/STO?® interfaces, 
thus opening the way to voltage-driven superconducting 
devices; it exhibits magnetic propertiesP® p, displays a 
strong and tunable 11 12 Rashba spin-orbit coupling? 3 it 
is extremely two-dimensional, having a lateral extension 
« 5 nm, thereby enhancing the effects of disorder due to 
extrinsic and/or intrinsic 1 ^ sources. 

Magneto-transport experiments reveal the presence of 
two kinds of carriers in LTO/STO, with high and low 
mobility, and superconductivity seems definitely to de¬ 
velop as soon as high-mobility carriers appear? 14 when 
the carrier density is tuned above a threshold value, 
by means of gate voltage V g . When the temperature 
T is lowered, the electrical resistance is reduced, and 
signatures of a superconducting fraction are seen well 
above the temperature at which the global zero resis¬ 
tance state is reached (if ever). The superconducting 
fraction decreases with decreasing V g , although a super¬ 
conducting fraction survives at values of V g such that 
the resistance stays finite down to lowest measured tem¬ 
peratures. When V g is further reduced, the supercon¬ 
ducting fraction eventually disappears, and the 2DEG 
stays metallic at all temperatures, seemingly undergo¬ 
ing weak localization at low T. At yet smaller car- 
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FIG. 1. Scheme of an oxide interface (not in scale), in the back 
gating configuration. The LAO (or LTO) topmost slab has 
a thickness of few nm, whereas the STO slab is « 0.5 mm 
thick. The back gate voltage V g is employed to tune the 
carrier density of the 2DEG. The top gating configuration 
is also possible. 


rier densities, the system behaves as an insulator. The 
width of the superconducting transition is anomalously 
large and cannot be accounted for by reasonable su¬ 
perconducting fluctuations?- 5 This phenomenology sug¬ 
gests instead that an inhomogeneous 2DEG is formed 
at these oxide interfaces, consisting of superconducting 
“puddles” embedded in a (weakly localizing) metallic 
background, opening the way to a percolative supercon¬ 
ducting transition. 16 Inhomo gene ities are revealed in var¬ 
ious magnetic experiments,® 710 in tunneling spectra?^ 
and in piezoforce microscopy measurementsP-® Seemingly, 
inhomogeneities at the nanometric scale coexist with 
larger (e.g., micrometric) scale inhomogeneities, revealed 
by the occurrence of striped textures in the current 
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distribution®^ and in the surface potential.^ 

Various aspects of the phenomenology of oxide inter¬ 
faces (henceforth referred to as LXO/STO interfaces, 
when referring to both LAO/STO and LTO/STO) have 
been separately discussed before. ® 1 ^ 5 h 6 E 1122I Here, we 
first provide new compelling evidence of the inhomoge¬ 
neous character of the 2DEG, extending previous multi¬ 
carrier analyses of magneto-transport measurements to 
deal with inhomogeneous systems. We then put together 
the various pieces of the jigsaw into an overall coherent 
theoretical framework. 

The plan of the present piece of work is the following. 
In Sec.[TlJ we propose a model for multi-carrier magneto¬ 
transport in inhomogeneous systems, and show that pre¬ 
vious analyses of the magnetoresistance and Hall resis¬ 
tance measurements in terms of two different species of 
carrier^® is fully consistent with the inhomogeneous char¬ 
acter of the 2DEG at LAO/STO and LTO/STO inter¬ 
faces. In Secs. |III[ |IV[ and [Vj we revisit some of our 
previous results. Assuming inhomogeneity as an em¬ 
pirical evidence, we show that resistance measurement^® 
and the topographic mapping of the superfluid density 10 
can be accounted for within a percolative scheme. In 
Sec. VI we discuss some preliminary aspects of a the¬ 
ory for metal-insulator-superconductor tunneling in in¬ 
homogeneous superconductors that is apt to r eproduce 
the measured tunneling spectral In Sec. VII we pro¬ 
vide further evidence for inhomogeneous superconductiv¬ 
ity at oxide interfaces, coming from the peculiar multiple 
quantum critical scaling, observed when superconductiv¬ 
ity is suppressed by means of a magnetic field, revisiting 
the results of Ref. [2T]- In Sec. VIII relying on the re¬ 
sults of Ref. [22], we show that the properties of the 
superconducting puddles (e.g., their fraction, and crit¬ 
ical temperatures) can be extracted from experiments 
and used to model intra-puddle multi-carrier supercon¬ 
ductivity, gaining insight about the pairing mechanism. 
Although some features of the diamagnetic response are 
seemingly related to strong superconducting coupling^® 
we show that inhomogeneities and multi-carrier super¬ 
conductivity fully account for the behavior of these sys¬ 
tems within a standard weak coupling BCS scheme. Con¬ 
cluding remarks are found in Sec.|IX| 


II. MULTI-CARRIER MAGNETO-TRANSPORT 
IN INHOMOGENEOUS SYSTEMS 

The detection of two species of carriers, with high and 
low mobility, by means of magneto-transport measure¬ 
ments in LTO/STO,® is not necessarily a direct evidence 
of inhomogeneity, since the two species could coexist in a 
homogeneous manner. However, to account for the phe¬ 
nomenology of the superconducting state of the 2DEG at 
LXO/STO interface, we proposed that the 2DEG is in¬ 
homogeneous, with higher density regions (the supercon¬ 
ducting puddle s) and lower density regions (the metal¬ 
lic background) .121123 Assigning a band structure to the 



FIG. 2. Hall resistance as a function of magnetic field for dif¬ 
ferent gate voltages V g , measured at T = 4.2 K in a LTO/STO 
15-unit-cell thick sample (sample A of Ref. [I]). Solid lines 
correspond to experimental data, taken from Ref. [3j, while 
the dashed lines are obtained here, fitting the data with Eq. 
@. For clarity, the curves have been shifted upwards by 
50 Q in increasing voltage order. Inset: Sheet resistance R, 
obtained via Eq. (JTJ , as a function of the gate voltage. 


system, one is then led to assume that, when the den¬ 
sity is large enough to fill the bands occupied by the 
high-mobility carriers, these should locally coexist with 
low-mobility carriers occupying the lower-lying levels (see 
also Sec. VIIII. A picture then emerges in which the low- 
mobility carriers alone are present in the metallic back¬ 
ground, whereas they coexist with the high-mobility car¬ 
riers in the superconducting puddles. To improve our 
description of LXO/STO interfaces, we must then rely 
on a theory for multi-carrier magneto-transport that is 
apt to deal with an inhomogeneous system. One such 
theory has been developed in the form of an Effective 
Medium Theory (EMT) for the Hall conductance of a bi¬ 
nary medium resulting from the mixture of two phases, 
based on rotation transformations (see Ref. [21] for a de¬ 
tailed description of this method). 


In our description of the LXO/STO interfaces, one 
(less dense) phase hosts one species of carriers, with 
low mobility, and the other (denser) phase hosts two 
species of carriers, with low and high mobility. Following 
Ref. [23], we first define the conductivity tensor for each 
of the two coexisting phases, in the presence of a mag¬ 
netic field of amplitude B. The diagonal elements of the 
conductivity tensor are 


aW(B) = _^_ 

xxK f l + [/3i(H)] 2 ’ 
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FIG. 3. (a) Fraction w of the high-mobility metallic phase 
[phase (2), see text]; (b) Densities of the less (m, empty cir¬ 
cles) and more (ri 2 , empty squares) mobile carriers extracted 
from the fits of the Hall resistance, Eq. (|2|, as a function of 
gate voltage measured at T = 4.2 K. (c) Mobilities /o (empty 
circles) and /12 (empty squares) of the majority (low-mobility) 
and minority (high-mobility) carriers, respectively, extracted 
from Eq. Type 1 (low-mobility) carriers only are present 
in the phase (1), while both type 1 and type 2 (high-mobility) 
carriers are present in the phase (2). 


Hereafter, the superscripts (1) and (2) label the two 
phases, the subscripts 1 and 2 label the two species of 
carriers (with low and high mobility, respectively), and 
we adopt the notations (3i(B) = /liB/c and cq = enj/Xj, 
where rii is the carrier density and /ii is the mobility of 
the i-th species of carriers. We point out that the peculiar 
aspect of our present description is that both species of 
carriers contribute in parallel to the conductivity tensor 
in the phase (2), which we identify with the higher den¬ 
sity superconducting puddles. We indicate with w = e 
the fraction of the system occupied by the superconduct¬ 
ing puddles, e being the deviation from the percolation 
threshold, that changes sign when the minority phase 
(w < 2 , e > 0) percolates and becomes the majority 
phase (w > e < 0), the threshold being w = \ in 
two-dimensional systems. 

Within EMT, the diagonal element of the conductivity 
tensor of the inhomogeneous (two-phase) system can be 
found, that correctly reproduces the limiting case of a 
pure phase, and takes the fornPH 
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Then, exploiting duality relations that connect the var¬ 
ious elements of the conductivity tensor when the minor¬ 
ity and majority phases are interchanged (see Ref.Bijl. 
one obtains the EMT expression for the off-diagonal el¬ 
ement of the conductivity tensor of the inhomogeneous 
(two-phase) system 
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From the above relations, the expression for the Hall re¬ 
sistance of the inhomogeneous system can finally be de¬ 
rived, 


r EMT 


i?f MT (S,e) = 


(£,< 


[ <J xy VlT ( B > e )] 2 + [ a xx^ T ( B ^ e )] 2 ' 


( 2 ) 


which we used to accurately fit the experimental Hall re¬ 
sistivity data of LTO/STO,® under strong magnetic field 
and at different gate voltages (see Fig. [ 2 ]). This proce¬ 
dure allows to extract the values of the mobilities ^ 
and of the densities n.j of the two species of carriers, as 
well as the fraction of the system occupied by the su¬ 


perconducting puddles, w. These values are reported in 
Fig.0 while the inset of Fig. [2] displays the evolution of 
the normal-state sheet resistance, extracted via Eq. <[T|) 
as R = [&xx IT ( b = 0,e)] -1 , with changing V g . 

Our present findings cast in a somewhat different per¬ 
spective the analysis carried out in Ref. [4] , where the 
appearance of more mobile carriers around V g = 0 V was 
found. While on the one hand we fully confirm that re¬ 
sult, on the other hand we find here that the magneto¬ 
transport data are well described assuming that these 
carriers do not appear uniformly in the whole system. 
Rather high- and low-mobility carriers coexist in a high- 
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FIG. 4. (a) Sketch of the behavior of the resistance as a function of temperature in a case when the superconducting fraction 
does not percolate, (b) Distribution of critical temperatures W(T C ) in the superconducting puddles, that occupy a fraction 
w of the sample (the remaining 1 — w fraction will never become superconducting), (c—g) When the temperature is reduced 
[from (g) to (c)] superconducting puddles appear in the system as soon as the local critical temperature exceeds T. However, 
if the puddles do not percolate down to T = 0, the global zero resistance state is never reached, (h) Sketch of the tunneling 
spectra measured in LAO/STO, as a function of bias voltage V : At temperatures well above the temperature at which the zero 
resistance state is reached (if ever), a suppression is observed in the tunneling spectra around zero bias, accompanied by more 
or less pronounced coherence peaks above the gap. The coherence peaks are smeared and the intra-gap spectra are broadened, 
with non-vanishing zero bias spectral weight, that we attribute to the metallic background. The curves are labelled by colored 
triangles that refer to the corresponding temperatures in panel (a). 


mobility metallic phase, here identified with the phase 
( 2 ), which is inhomogeneously distributed at the interface 
and spatially separated from the less metallic phase ( 1 ), 
where only one species of less mobile carriers is present. 
In this framework, the enhanced conductivity around 
V g = 0 V and the changes of slope of the Hall resistance 
at high magnetic field occur because a finite fraction w of 
phase ( 2 ), characterized by an overall higher carrier den¬ 
sity and hosting high-mobility carriers, appears around 
this gate voltage [see Fig.[3](a)]. In agreement with the 
conclusions of Ref. [J] , we also find that the more mobile 
carriers have a lower density than the low-mobility ones 
[Fig.[3](b)] but, at odds with the results obtained within a 
homogeneous description, in the present inhomogeneous 
scheme the mobilities of the two species stay rather con¬ 
stant over the whole range of gate voltages [see Fig.|3jc)] . 


III. PERCOLATIVE SUPERCONDUCTIVITY 

After assessing the occurrence of inhomogeneous multi¬ 
carrier magneto-transport at LXO/STO interfaces, we 
discuss the superconducting transition that is driven by 
tuning the gate voltage V g (i.e., the carrier density) above 
a threshold value. It has been shown that the supercon¬ 
ducting transition occurring in inhomogeneous systems 
is well described within the EMT.0I3 EMT is a mean- 
field-like theory apt to describe a random resistor net¬ 
work (RRN) that lacks spatial correlations. The EMT 
equations are obtained embedding one given random re¬ 
sistance Ri in an effective medium of constant resistance 
R. This latter is chosen in such a way as to have the same 
current flowing through Ri as in the RRN. The EMT re¬ 
sistance can be showrP^l to be larger than the parallel 
and smaller than the series of the random resistances, 


the two limiting values being reached in infinite and one 
dimensions, respectively. 

The resistance of the LTO/STO interface exhibits a 
marked suppression due to incipient superconductivity, 
that is accurately fitted^ assuming that the supercon¬ 
ducting puddles occupy a fraction w < 1 of the sample, 
and that each puddle is characterized by a random local 
critical temperature T c . For the sake of definiteness, we 
adopt a Gaussian distribution of critical temperatures, 
W ( T c ), parametrized by its mean value T c and its width 
7 . The remaining 1 — w fraction of the sample is occupied 
by the metallic background. The resistan ce at tempera¬ 
ture T is found within EMT to be 15 lb 22 


R(T) = R 0 


(1 — w) + w erf 


T - T, 
lV2 


and results form the metallic background (first term in¬ 
side the square brackets) and from not yet superconduct¬ 
ing puddles (i.e., those puddles with T c < T, second term 
inside the square brackets, erf being the error function); 
the remaining puddles (i.e., those with T c > T), have 
become superconducting and do not contribute to the 
resistance. The high-temperature resistance i?oo, w, T c , 
and 7 are used as fitting parameters. The global zero 
resistance state is reached at the percolative transition 
temperature T p < T c such that 


erf 


( Ty-Tc \ 

\ lV2 J 


w — 1 
w 


(3) 


A solution for T p only exists if the superconducting frac¬ 
tion of the 2DEG can percolate in the two-dimensional 
system, i.e., for w > When T p < 0, or when it is not 
at all defined (for w < |), the resistance remains finite 
down to T = 0, although the presence of a sizable (yet 
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FIG. 5. (a) Measured sheet resistance in a LAO/STO sample, 
for V g = -10,10, 30, 50, 70,100,150, 200 V. The solid lines are 
the EMT fits, (b) Same as in panel (a), with fewer experi¬ 
mental curves, fitted by RRN curves (solid lines), (c) Weight 
of the superconducting fraction and (d) width of the T c dis¬ 
tribution, from EMT (line and small filled squares) and RRN 
(large open squares) fits. 


not percolating) superconducting fraction is mirrored by 
a sizable suppression of R(T), as sketched in Fig.[4](a). 

The resistances measured in a LTO/STO sample as 
a function of T for various values of the gate voltages 
V g Pare shown in Fig.^a) together with the fitting EMT 
curves (symbols and solid lines, respectively). Fig.[5]is in¬ 
spired to a similar figure of Ref. [22] , with new elements 
included to make explicit contact with the overview con¬ 
tained in this piece of work. The set of parameters 
Rao,w, T c , 7 , i.e., the distribution of the puddles, changes 
with changing the carrier density by means of gating. 

We are thus able to extract from the fits the mean 
intra-puddle critical temperature T c [i.e., the temper¬ 
ature at which R(T) changes curvature within EMT], 
which will be analyzed in Sec. |VIlH and the overall super¬ 
conducting fraction w, tracked by the solid line with filled 
squares in Fig.[5jc) . We point out that EMT disregards 
spatial correlations, so that the presence of pronounced 
tails in the resistance, in the presence of a percolating 
superconducting cluster, force the overall superconduct¬ 
ing fraction to bc 1 ^ to ft |, as shown in Fig. pic). The 
width 7 of the Gaussian distribution of T c is plotted as a 
solid line with filled squares Fig.[5jd) . It increases as the 
fraction w of the puddles goes to zero. This is rather nat¬ 
ural because a reduction of the carrier density is expected 
to emphasize the effects of disorder, so that fluctuations 
of the local superconducting critical temperature should 
increase, leading to a broadening of the T c distribution. 

A comment is now in order. Dealing with supercon¬ 
ducting puddles embedded in a metallic background, we 
expect the proximity effect! 2 ^ to play an important role in 
LXO/STO interfaces. Within EMT, this effect sure en¬ 
tails a temperature dependence of the superconducting 


fraction w. However, when fitting the resistance curves 
R(T), the introduction of a temperature dependent w(T) 
is not viable, because it would made a good fit by def¬ 
inition. Nonetheless, the role of proximity effect can be 
analyzed when discussing tunneling spectra. We shall 
come back to this point in Sec. [Vi] 


IV. SUPERFLUID DENSITY IN 
INHOMOGENEOUS SUPERCONDUCTORS 


The inhomogeneous character of the 2DEG at the 
LXO/STO interfaces raises the question of the descrip¬ 
tion of the superfluid properties in a mixture of two 
phases. The superfluid density n s was measured 10 ^ by 
means of a SQUID in LAO/STO interfaces. The mea¬ 
surements average over micrometric scales and are there¬ 
fore not sensitive to submicrometric inhomogeneities. 
Nonetheless, the idea that the 2DEG at these interfaces 
is inhomogeneous, is still supported by the evidence of 
variations of n s on the micrometric scale within a given 
sample. Encouraged by the marked similarity of the re¬ 
sistance curves in LAO/STO and in LTO/STO, we apply 
EMT also to describe the measured local n s as an average 
over an inhomogeneous state of submicrometric puddles. 
We point out that recent experiments in LTO/STC^Sled 
to estimate the typical size of the puddles to be ft 100 ran, 
thereby supporting this idea. 

We proceeded extending EMT to small frequency 
and assigning to the metallic background a Drude-like 
complex conductivity <7m{w) = A (r + iw ) -1 and to 
the superconducting puddles a purely reactive conduc¬ 
tivity crs(w) = A(lw) -1 . We then defined the resistivity 
AsM = cr^w) and p M {w) = = p 0 + AsM, 

with po = r/A. At high T, the system is metallic and 
p(ui) = pm(w). However, when the temperature is low¬ 
ered, the static resistivity vanishes within each individual 
puddle as soon as T equals the local T c . Although the full 
expression of the complex resistivity can be found within 
EMT, aiming at describing the static diamagnetic re¬ 
sponse, we give only the expression up to terms ~ ui, i.e., 
p(w) « po (w M - w s ) w M - w s ) + | w M - w s \~ 1 p s (u ;), 
where A is the Heaviside function, ws is the fraction of 
puddles that have become superconducting (at a given 
temperature), and wm = 1 — ws is the non superconduct¬ 
ing fraction (resulting both from puddles that have not 
yet become superconducting and from the metallic back¬ 
ground). Evidently, when wm > ws, the conductivity is 
Drude-like. However, below the percolation temperature 
T p (whenever defined), wm < ws, and the conductivity 
is purely reactive cr(ui) = T (ws — wm) (*w) _1 . There¬ 
fore, employing also Eq. ^ , we find that the superfluid 
density of the percolating two-dimensional network for 
T < T p is 


n s oc ws — w M 


erf 


T p -T c 

7^2 


— erf 


T-T c 

7v/2 
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FIG. 6. (a) Superfluid density (symbols) as a function of T 
and EMT fits according to Eq. Q (solid lines). Weight of the 
superconducting fraction (b) and width of the T c distribution 
(c) extracted from the EMT fits, as a function of V g . The 
lines are guides to the eye. 


of 7 resembles that obtained from transport in LTO/STO 
[Fig.[5](d)] and is of comparable magnitude. 

Of course, a proof of the above arguments requires a 
model accounting for space correlations, that includes 
both closed loops (relevant for diamagnetism) and con¬ 
nected paths (relevant for transport). While such a 
model has been found for transport (see Sec.[v]), a simi¬ 
lar model for the diamagnetic response is not yet avail¬ 
able, so our discussion on this specific aspect has a purely 
speculative character. An experimental test to our sce¬ 
nario could be provided by the observation of differ¬ 
ent diamagnetic responses in field-cooled and zero-field- 
cooled samples. In the first case one expects substantially 
smaller diamagnetic fraction, since a sizable flux would 
be trapped in the normal part encircled by the super¬ 
conducting loops (see the shaded region in Fig. [f]and the 
discussion in Sec.[V|. 

The study of the dynamical response of an inhomoge¬ 
neous superconductor at finite frequency, within EMT, 
is presently under investigation. A somewhat similar 
approach, applied to superconducting stripes in high- 
temperature superconductors, can be found in Ref. |28l . 


In Fig.pl a), this average superfluid density and our 
EMT fits [Eq. 0] are reported, for different values of V g . 
Fig .[ 6 ] is inspirea to a similar figure of Ref. [22], again suit¬ 
ably modified to make explicit contact with the present 
discussion. The behavior of n s (T), although qualitatively 
resembling the BCS prediction, may quantitatively differ 
from it. The slope at T = T p , for instance, is controlled 
by the width 7 of the distribution W(T C ). Thus, the de¬ 
viations from standard BCS prediction measured in Ref. 
;IOj. and attributed to a tendency to strong coupling, 
are here alternatively explained claiming that the 2DEG 
is in a weak coupling BCS regime (see Sec. VIII), but 
the superfluid density is intrinsically inhomogeneous at 
the submicrometric scale and is averaged at the micron 
scale by the SQUID pick-up loop used in Ref. [ID] . We 
point out that our mean field approach is justified by 
the observation that the channels down to 500 nrn wide 
share properties similar to those of the larger samples.^ 
Our fits yield the fraction w of volume occupied by the 
puddles, extracted from measurements of the superfluid 
density and reported in Fig. Glb). It ranges from | to 
1 , and is always larger than me fraction obtained from 
transport measurements. This is not surprising, because 
transport mainly probes the long-range connectivity of 
the percolating path, regardless of dead ends and discon¬ 
nected superconducting regions,^, whereas diamagnetic 
screening measurements are sensitive to all sufficiently 
large superconducting loops, even when not connected to 
the backbone. In this case, the diamagnetic fraction can 
be large, while the long-distance connectivity is small, 
if many puddles or loops are disconnected (see Sec. 0 - 
Fig-i c) reports the behavior of the width 7 of the distri¬ 
bution W(T C ), inferred from the superfluid density. We 
point out that, despite the fact that we are dealing with 
different materials and physical quantities, the behavior 


V. SPACE CORRELATIONS WITHIN THE 
SUPERCONDUCTING CLUSTER 

A possible explanation for the discrepancy between the 
superconducting fraction w observed in transport and 
diamagnetism measurements may rest upon the filamen¬ 
tary structure of the superconducting cluster at oxide 
interfaces.^] So far, we have made use of the mean-field¬ 
like EMT, which completely neglects spatial correlations. 
To investigate the mechanisms determining the supercon¬ 
ducting fraction observed in transport, w « 1 , we solved 
a RRN where the superconducting puddles form a spa¬ 
tially correlated cluster embedded in a metallic matrix. 
Preliminary results 1 ® indicated that a superconducting 
cluster which is dense at short distances and filamentary 
at larger distances is necessary in order to reproduce the 
observed tails of the resistance curves near percolation. 

To this purpose, we generated a fractal-like cluster 
with small long-scale connectivity, that percolates only 
when almost all bonds have become superconducting. We 
point out that the fractality of the clusters is an artifice 
to produce spatially correlated networks that are at the 
same time dense at short distances and filamentary over 
long distances. For a Gaussian W(T C ), the low tempera¬ 
ture tail of the distribution must be necessarily reached, 
and a correspondingly pronounced tail in the resistance 
is obtained. Compact clusters fail to reproduce the tails 
in the resistance, in the absence of filamentary struc¬ 
tures over long distances. Our systematic investigation 
showed that the presence of loosely connected filaments is 
a necessary ingredient to reproduce the behavior of R(T) 
near percolation. On the other hand, a purely filamen¬ 
tary structure, no matter how dense at short distance, 
is too loose and is not apt to describe the behavior of 











VI. TUNNELING SPECTRA OF 
INHOMOGENEOUS SUPERCONDUCTORS 
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FIG. 7. Sketch of the RRN with superpuddles (green regions) 
connected by filaments of puddles (purple circles). The large 
shaded region, that does not contribute to superconducting 
transport, contributes instead to the diamagnetic response in 
the case of a zero-field-cooling measurement, since the super¬ 
conducting currents confined in the superpuddles and loop 
structures encircled by the puddles are apt to screen a mag¬ 
netic field. In the case of a field-cooling measurement, one 
expects instead that the magnetic flux is frozen within the 
loop structures encircled by the puddles. 


R(T) at higher temperatures. In order to tune the den¬ 
sity of the superconducting cluster without significantly 
changing the long-distance connectivity, we decorate the 
filaments with randomly distributed superpuddles, for 
simplicity assumed circular, their number and size be¬ 
ing chosen to produce weights w ranging from 0.3 to 0.7 
(see the sketch in Fig.[7|. Superpuddles may be produced 
by by large tails in the distribution of the puddle sizes, 
or by extrinsic pinning centers, promoting the nucleation 
of much larger puddles. We systematically investigated 
the effect of the size and density of the superpuddles. A 
smaller fraction of larger superpuddles or a larger fraction 
of smaller superpuddles are essentially equivalent, as long 
as the superpuddles do not overlap to form percolating 
clusters. 

In Fig.[§b), we show our fits of the resistance obtained 
within our RRNs. Noticeably, even though the super¬ 
conducting fraction w is no longer forced to |, as it was 
instead within EMT, the resistance displays pronounced 
tails only if 0.50 < w < 0.65. The lower bound is imposed 
by the high slope at intermediate temperatures whereas 
the upper bound is due to the pronounced tail near per¬ 
colation. In Figs.[5jc) and (d) we show the fraction w 
occupied by the superconducting cluster and the width 
7 of the Gaussian distribution W ( T c ) obtained within 
the RRN (open circles): 7 qualitatively resembles the 
width obtained within EMT, and increases upon lower¬ 
ing V g . Thus, EMT and RRN models lead to similar re¬ 
sults about the distribution of T c and its variation with 
V g , provided the RRN is dense at short distances and 
filamentary over long distances. 


Another possibility of probing the inhomogeneous 
character of the supercondutcting state at LXO/STO 
interfaces is provided by tunneling experiments. Re¬ 
cent metal-insulator-superconductor tunnel spectroscopy 
measurements on LAO/STOj^reveal the occurrence of a 
state with finite resistance, but superconducting-like den¬ 
sity of states (DOS). The measurements are performed 
depositing a metallic Au electrode on top of the insulat¬ 
ing LAO layer, and applying a bias voltage V to drive a 
tunnel current I between the electrode and the 2DEG. 
The electrode measures several hundreds /im across, thus 
being orders of magnitude larger than the nanoscale 
inhomogeneities.^ At the lowest measured temperature, 
T = 30 mK, the spectra reveal a gap in the DOS at the 
Fermi energy, over the entire range of explored gate volt¬ 
ages V g € [—300,300] V, accompanied by more or less 
broadened coherence peaks above the gap, pointing to 
superconducting coherence and pairing as the origin of 
DOS suppression [see Fig.Qh)]. In the carrier depleted 
regime ( V g <C 0 ), the suppression is found even in the 
absence of global superconductivity, again highlighting 
the inhomogeneous character of the state formed by su¬ 
perconducting puddles embedded in a metallic matrix. 
At very low carrier concentration, V g ~ —300 V, the co¬ 
herence peaks have been completely smeared, although a 
sizable gap is still present as a signature of (incoherent) 
pairing. We shall call pseudo-gap state a state display¬ 
ing at the same time a DOS suppression in the tunneling 
spectra, accompanied by more or less pronounced coher¬ 
ence peaks, and non-zero electrical resistance. At higher 
carrier concentration, V g = 200 V, the 2DEG displays a 
superconducting gap and coherence peaks that decrease 
with increasing temperature and vanish around 300 mK 
(which agrees with the critical temperature of bulk STO 
reported in Ref. [22])• At V g < 0, the gap closes and 
the coherence peaks vanish, at temperatures much higher 
than the global T c . 

Since, as it was observed above, the tunneling spec¬ 
tra are taken over a spot measuring several hundreds /im 
across, we proposed^ that the observed pseudo-gap re¬ 
sults from an average over superconducting regions (the 
superconducting puddles, with a DOS described by stan¬ 
dard BCS theory) and metallic regions (composed of the 
metallic background and of not yet superconducting pud¬ 
dles) with constant DOS N 0 . 

The differential conductance is usually written as 

d T f°° 

— (V)=G 0 ~G 1 J f'(E + eV)N(E)dE, (5) 

where N(E) is the DOS, f'(E) is the derivative of the 
Fermi distribution function, the positive constant Go cus¬ 
tomarily accounts for effects such as leakage currents, and 
G\ is a dimensional constant. We recall that in our model 
superconducting pairing occurs within each puddle below 



a local critical temperature T c , randomly distributed ac¬ 
cording to a probability distribution W(T c ). The DOS of 
the 2DEG probed in tunneling spectra has three distinct 
contributions, 

N(E) = (1 - w)N 0 + wN 0 f d T c W{T C ) 

J — OO 

/ oo 

d T c W(T c )N a( t c ,t)(E). ( 6 ) 

The first two terms correspond to the metallic back¬ 
ground and to puddles where pairing has not taken place 
yet, respectively, and fully account for the zero bias back¬ 
ground observed in the tunneling spectra, allowing us to 
take G 0 = 0 in Eq. <©• The third term corresponds to 
puddles that developed a finite pairing gap A. The DOS 
within these puddles is taken as 


N a (E) = 


( 1 -x) 


\E | 


y/E 2 - A 2 


N o 0(\E\- A). 


The first term is the standard BCS expression, and de¬ 
scribes coherent pairing occurring within a (1 — x) frac¬ 
tion of the gapped part. The second term describes 
puddles that, although having a gap, are too small to 
exhibit phase coherence. We define w pa i r as the total 
fraction of the system that can display pairing down to 
T = 0, and the coherently and incoherently paired frac¬ 
tions w co h = (1 — x) Wp a i r and Wi nc = x w pa i r . The latter 
term is necessary because the experimental spectra are 
gapped but display no coherence peaks when V g <C OPS 
The gap is assumed to take the BCS expression 

A (T C ,T) = 1.76 T c tanh 

The value of No is readily determined by the high-bias 
part of the spectra, while w , x. and the parameters of 
the distribution of critical temperatures, T c , 7 , are used 
as fitting parameters. 

At low temperature, accurate fits the spectra and of 
their evolution as a function of the gate voltage V g were 
obtained.^ Remarkably, the width of the distribution of 
critical temperatures, 7 , turned out to be very weakly 
dependent of V g , indicating that the distribution of T c in 
the sample is essentially related to structural properties, 
like the local amount of disorder. 

When fitting the temperature dependence of the spec¬ 
tra at fixed V g , the attempt to make use of the same 
temperature-independent set of parameters, although 
capturing the main features of the spectra, yielded 
fits definitely less convincing than the fits at low tem¬ 
perature. However, releasing the severe constraint of 
temperature-independent parameters, very good fits are 
obtained by letting w and x vary, while the parameters of 
the distribution of critical temperatures, T c , 7 are kept 
fixed. The variation of w and x with the temperature 
is more clearly understood when expressed through the 
quantities w co h and w co h defined above. It turns out that 


the fraction of the sample occupied by the superconduct¬ 
ing cluster increases with decreasing temperature, and 
saturates at low T, likely indicating that a sizable part 
of the metallic background is g radually driven supercon¬ 
ducting by proximity effect 


VII. MULTIPLE QUANTUM CRITICAL 
BEHAVIOR IN INHOMOGENEOUS 
SUPERCONDUCTORS 


The recent analysis of the quantum critical behavior 
that is observed when the superconducting phase is sup¬ 
pressed by means of a magnetic field perpendicular to 
the interface ;- 1 provided further support in favor of our 
scenario for inhomogeneous superconductivity at oxide 
interfaces. When the resistance R of LTO/STO, rescaled 
by a characteristic value R, is plotted as a function of 
the variable ( B — B)/T 1 ^ zv , where B is a characteristic 
magnetic field (see Fig. [ 8 ]), z is the dynamical exponent 
(converting a critical length scale in a time scale), and v 
is the critical exponent of the correlation length as a func¬ 
tion of temperature, two quantum critical scaling regimes 
are found, in different temperature ranges. These two 
scaling regimes, separated by a crossover at intermediate 
temperatures, are in correspondence with two different 
values of the characteristic magnetic field B = B X ,B C . 
The smaller field, B x , is related to the scaling at higher 
T, with zv ss |, whereas the (slightly) larger field, B Cl 
corresponds to the scaling at lower T, with zv k, |. No¬ 
ticeably, the characteristic fields B x and B c coalesce at 
low carrier density (i.e., low gate voltage V g ). When V g is 
increased, B x saturates to a constant, whereas B c closely 
tracks the superconducting critical temperature T c (with 
a conversion factor of 1 T corresponding to IK). 

A possible explanation for this multiple critical be¬ 
havior relies on the assumption that superconductiv¬ 
ity within an isolated puddle would be suppressed by 
the smaller critical field B x . However, the puddles 
are coupled, being embedded in a common metallic 
background . 31 When inter-puddle coupling eventually in¬ 
tervenes, superconductivity is strengthened and survives 
up to the (slightly) larger critical field B c . If we bor¬ 
row the value of the dynamical exponent z = 1 , which 
is commonly adopted in similar situations,®^ the intra¬ 
puddle criticality is described by the XY modeP- in the 
clean limit, where v ss |. However, with lowering the 
temperature, the coherence length increases, eventually 
exceeding the puddle size, and inter-puddle superconduc¬ 
tivity establishes in the inhomogeneous landscape of pud¬ 
dles embedded in the metallic background. In this case, 
the exponent v must obey the Harris criterion 34 for dis¬ 
ordered systems, and indeed we find v k, | > 1 . 








9 




FIG. 8 . Sketch of the procedure adopted to extract the quan¬ 
tum critical behavior when the superconducting phase is sup¬ 
pressed by means of a magnetic field perpendicular to the 
interface, (a) When the resistance R is plotted as a function 
of temperature T, three behaviors are observed over a certain 
temperature interval, exemplified by the curves corresponding 
to the three magnetic fields B\ < B 2 < B 3 : superconducting, 
critical, insulating, (b) To better identify the critical value 
of the magnetic filed, B x or B c , isotherms are plotted as a 
function of B, exemplified by the three curves corresponding 
to the temperatures Ti < T 2 < T 3 : the crossing point of the 
isotherms corresponds to B ( B x or B c ), and the correspond¬ 
ing characteristic value R is obtained on the resistance axis. 


VIII. MULTI-CARRIER BCS MODEL 


According to the discussion developed so far, the in¬ 
homogeneous character of the 2DEG at the LXO/STO 
oxide interfaces induces a distribution W(T C ) of local crit¬ 
ical temperatures, its mean value T c depending on the 
overall carrier density (i.e., on V g ). This dependence is 
obtained fitting the resistance data in Fig.[sja) within 
EMT, and sheds light on the intra-puddle pairing mecha¬ 
nism. As discussed in Sec.[TJ detailed magneto-transport 
measurements^ highlighted the coexistence of a sizable 
amount of low-mobility carriers and a smaller amount of 
high-mobility carriers in LTO/STO. Here, as well as in 
LAO/STOpSl superconductivity seems definitely to de¬ 
velop as soon as the high-mobility carriers appear. 

Accordingly, we proposed 22 that superconducting pair¬ 
ing within the 2DEG formed at the oxide interface may 
be described by a multi-bancP^®^ BCS-like Hamiltonian 

HbCS = ( a k,£,f a k,£,t + 4,<,i a k,<,l) 

M 

, ' 9w t 

+ jy a M,t a —M,4. a —k',f',.f a k'/',t D' 

k,£ 

k',£' 

where e CT (a k 1 a) creates (annihilates) an electron with 
two-dimensional wave vector k = ( k x , k y ), parallel to the 
plane of the interface, and spin projection er =f, j., be¬ 
longing to the £-th sub-band, with dispersion law 


£k,f 


Vkl | hW y 


where m^ x < y \ are the (possibly anisotropic) effective 
masses of the charge carriers and g is the chemical po¬ 
tential, gw are the intraband (for t = £') or interband 



FIG. 9. Sketch of the quantized levels Ee in the potential well 
V(z) confining the 2DEG at the interface (z is the coordi¬ 
nate perpendicular to the interface), and of the corresponding 
carrier density distributions. The lower-lying levels, e.g., at 
energies E = £o,£i, accommodate the low-mobility carriers 
(LMC), whereas the topmost level, e.g., at E = £ 2 , hosts the 
high-mobility carriers (HMC). 


(for £ ^ £') pairing amplitudes and N is the number 
of k points within the first Brillouin zone. The sub- 
ban ds can originate, e.g., from the multiband structure of 
STO 12 41 and/or from the quantum confinement within 
the self-consistent potential well at the interface )^i2l i n 
our schematic description, we represent the whole set of 
low-lying bands with one sub-band (£ = 1) accommodat¬ 
ing the non superconducting low-mobility carriers, while 
the high-mobility carriers in the sub-band £ = 2 are 
paired and give rise to a finite T c (see Fig.[ 9 j where a 
sub-band structure originating from the quantum con¬ 
finement is sketched, readapted from a similar sketch in 
Ref. gU). 

Thus, according to our interpretation, the supercon¬ 
ducting puddles are regions where the £ = 2 sub-band 
is locally filled, whereas the (weakly localizing) metallic 
background corresponds to regions where the £ = 2 sub¬ 
band is empty. The phenomenology of the superconduct¬ 
ing phase at oxide interfaces is reproduced assuming that 
the pairing amplitudes are such that gn <C (312,321) "C 
322- This condition is also consistent with the analysis of 
a two-band model in Ref. [ 43 ;- For simplicity, to reduce 
to a minimum the number of free parameters, in the fol¬ 
lowing we assume that 311 = 312 = 321 = 0. According 
to the standard BCS approach, the pairing amplitudes 
are only effective in a window |£k,t|> ICk'.r I < kui 0, where 
ujq is a characteristic cut-off frequency. The prime su¬ 
perscript attached to the last sum in Eq. Q implies 
this restriction. We assume that the bottoms of the two 
sub-bands are well separated, £2 — £i Sw0, and take 
henceforth £2=0. 

In principle, when deducing a BCS-like Hamiltonian, 
one should take care of vertex corrections (which are 
expected to be relevant when Migdal’s condition is vio¬ 
lated). However, the task of deducing a BCS-like hamil- 
tonian in a multi-band model, where electrons with large 
and small Fermi energy coexist, is overwhelmingly diffi¬ 
cult. We rather take Eq. ([TJ) as a phenomenological low- 
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of work. The various band masses are taken all equal to 
0.7 electron masses, corresponding to a scenario where 
the sub-band structure originates from quantum confine¬ 
ment of the lowest t 2 g band of bulk STO near the inter¬ 
face. We then obtain the dimensionless coupling constant 
A ~ 0.125 (therefore, consistent with our assumption that 
the system falls in the weak coupling regime) and a cut¬ 
off energy scale Hloq ss 23 meV, which is compatible with 
a typical phonon energy in STOpUl 


IX. CONCLUDING REMARKS 


FIG. 10. Average critical temperature extracted from the 
EMT fit of the sheet resistance data (empty squares), and 
critical temperature of our BCS-like model (solid line with 
diamonds), as a function of V g . T c starts to rise when the 
chemical potential enters the HMC sub-band ( fj, = e 2 ) and 
saturates to the standard BCS value when q = e 2 + hu>o. The 
two regimes are separated by the dashed vertical line. 


energy effective Hamiltonian, imagining that the various 
high-energy effects have been accounted for by a suitable 
dressing of the bare physical parameters. Thus, the iden¬ 
tification of wo with the characteristic frequency of the 
pairing mediator is expected to hold only indicatively. 

For fj, < 0 the sub-band hosting the high-mobility car¬ 
riers is empty and T c = 0. For 0 < /i < fiwo, pairing 
occurs at the critical temperature 

T c ss 1.14 e -1//A , 

where A = <722 N BMC is the dimensionless superconduct¬ 
ing coupling and N BMC is the DOS of the sub-band filled 
by the high-mobility carriers. For n > hui 0 , T c saturates 
to the standard BCS value 

T c BCS « 1.14fiw 0 e“ 1/A . 

The previous results can be cast into a single expression 


T c (fj) = Tf cs d(/z) min 


hue) 


,1 


( 8 ) 


Remarkably, the amplitude A// of the interval in cor¬ 
respondence of which T c is an increasing function of /r 
provides a direct measure of the cut-off energy scale, 
Huiq = A/q which, with the caveat recalled above, we 
interpret as a possibly crude estimate of the characteris¬ 
tic energy scale of the pairing mediator. 

The fit of the curve T c (V g ), extracted from the ex¬ 
perimental data in Fig.[5](a) within EMT, with the curve 
T c (n) in Eq. ([8]), exploiting the approximately linear re¬ 
lation between V g and /q® yields the result illustrated in 
Fig. [lO] (orange line with diamonds and empty squares, 
respectively). Fig. 10 is inspired to a similar figure of 
Ref. j22j , redrawn here in a modified fashion, to make ex¬ 
plicit contact with the overview contained in this piece 


In conclusion, we first analyzed magneto-transport ex¬ 
periments in LTO/STO oxide interfaced within a multi¬ 
carrier EMT model, confirming the occurrence of two 
kind of carriers, with high and low mobility, within an 
inhomogeneous landscape. We proposed that the sys¬ 
tem consists of regions with higher carrier density, were 
both carrier coexist, and regions with lower carrier den¬ 
sity, where only the low-mobility carriers are present. We 
also confirm that the high-mobility carriers have a lower 
density than the low-mobility ones. The mobilities of the 
two species turn out to be almost independent of the gate 
voltage. Thus, the enhancement of conductivity observed 
around V g = 0 V and the change in the slope of the Hall 
resistance at high magnetic field occur because a finite 
fraction w of regions with higher carrier density, hosting 
high-mobility carriers, appears around this gate voltage. 

We then described superconductivity in LAO/STO 
and LTO/STO within a scenario in which superconduct¬ 
ing puddles (the regions with coexisting high and low- 
mobility carriers) are embedded in a metallic background 
(the regions with low-mobility carriers only), and form a 
percolating network. In this framework, the sheet resis¬ 
tance of LTO/STO interfaces is very well described by 
EMT or by a RRN for an inhomogeneous 2DEG with 
a substantial filamentary character. Fitting the exper¬ 
iments, we were able to extract the random distribu¬ 
tion of T c at at various V g (i.e., various carrier densi¬ 
ties ). A similar approach was adopted to fit the rnicro- 
metrically averaged superfluid stiffnesiP 0 1 and the pseudo¬ 
gap in tunneling spectrjPS of LAO/STO. In particular, 
we showed that our model accounts well for the metal- 
insulator-superconductor tunneling spectra measured in 
LAO/STO, and allows us to draw the conclusion that the 
fraction of the sample occupied by the superconducting 
puddles increases when the temperature is reduced, and 
saturates at low temperature, likely indicating that also 
a sizable fraction of the metallic background gradually 
becomes superconducting by proximity effect. 

Assuming an effective two-band model with supercon¬ 
ductivity triggered by the presence of few high-mobility 
carriers, locally filling the highest-energy band, we ac¬ 
count for the density dependence of the intra-puddle T c 
within a simple BCS weak coupling scheme. As an im¬ 
portant by-product, we find that the range of variation 
in V g of the (average) intra-puddle T c is directly related 
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(via the chemical potential /i) to the cut-off energy scale 
Huiq. Taking this value as a crude estimate of the typical 
energy scale of the pairing mediator, we find that this is 
compatible with phonon-mediated superconductivity. 

Although this was not the main focus of the present 
piece of work, few words are now in order to discuss the 
origin of electron inhomogeneity at LXO/STO interfaces. 
Extrinsic mechanisms® like impurities and growth de¬ 
fects are always accountable of rendering an interface in¬ 
homogeneous. However, experimental data show that in¬ 
homogeneity is never reduced below a still sizable extent, 
even in the best samples. Furthermore, the observation of 
negative electron compressibility in a low-filling regime— 
suggests that intrinsic mechanisms (e.g., in the form of 
effective electron-electron attractions) are present, which 
may render the 2DEG formed at the LXO/STO interface 
inhomogeneous by phase separation, even in a perfectly 
clean and expectedly homogeneous system. These facts 
led us to look for intrinsic mechanisms of inhomogene¬ 
ity. A possible intrinsic mechanism for the formation 
of an inhomogeneous 2DEG was proposed)- 121 ® that re¬ 
lies on the strong density dependent RSOC inferred at 
these interfaces)— and yields a non-rigid band structure, 
with the possibility that the chemical potential of the 
carriers is a non-monotonic function of the carrier den¬ 
sity for reasonable values of the model parameters, giv¬ 
ing rise to a thermodynamic instability. The resulting 


RSOC-driven electronic phase separation provides a nat¬ 
ural framework for the occurrence of regions with higher 
and lower carrier density (the superconducting puddles 
and the metallic background, respectively), accounting 
for the inhomogeneous character of the superconducting 
phase at oxide interfaces. Of course, whenever electronic 
phase separation comes about, one should worry about 
the Coulomb energy cost to be paid for the segregation of 
charge carriers. Thanks to the large dielectric constant of 
STO, we showed® that the size of the domains can grow 
as large as several tens of nanometers. An alternative 
route to intrinsic inhomogeneity might be provided by 
the pronounced density dependence of the self-consistent 
potential well V(z) confining the 2DEG at the interface 
(see Fig. [9]). Calculations within a Schrodinger-Poisson 
self-consistent scheme® show indeed that also this mech¬ 
anism is apt to yield a non-monotonic chemical potential, 
again resulting in electronic phase separation. 
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